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1 METHOD AND SYSTEM FOR COMPUTING ROAD GRADE DATA 

2 REFERENCE TO RELATED APPLICATION 

3 The present application is related to the co-pending application entitled 

4 "METHOD AND SYSTEM FOR OBTAINING ROAD GRADE DATA" filed on the 

5 same date herewith, Attorney Docket No. N0154US, the entire disclosure of which is 

6 incorporated by reference herein. 
7 

8 BACKGROUND OF THE INVENTION 



9 The present invention relates generally to geographic databases, and more 

10 particularly, the present invention relates to a method and system for computing road 

1 1 grade data using collected altitude or elevation data. 

12 Geographic databases have various uses. Geographic databases are used in 



13 in- vehicle navigation systems, personal computers, networked computing environments, 

14 and various other kinds of platforms, as well as on the Internet. Geographic databases are 

15 used with various kinds of applications to provide various navigation-related and map- 

16 related functions including map display, route calculation, route guidance, truck fleet 

17 deployment, traffic control, traffic monitoring, electronic yellow pages, roadside 



18 assistance, emergency services, and so on. 

19 In order to provide these kinds of functions, a geographic database includes data 

20 that represent geographic features in a region. The geographic features that are 

21 represented in a geographic database may include roads, intersections, and so on. A 

22 geographic database includes information about the represented geographic features, such 

23 as the geographic coordinates of roads in a geographic region, speed limits along the road 

24 segments, locations of stop lights, turn restrictions at intersections of roads, address 

25 ranges, street names, and so on. 

26 Collecting information for a geographic database is a significant task. Not only is 

27 the initial collection of data a significant undertaking, but a geographic database needs to 
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1 be updated on a regular basis. For example, new streets are constructed, street names 

2 change, traffic signals are installed, and turn restrictions are added to existing roads. 

3 Also, new levels of detail may be added about geographic features that are already 

4 represented in an existing geographic database. For example, an existing geographic 

5 database for roads may be enhanced with information about lane widths, shoulder sizes, 

6 traffic signs, lane barriers, address ranges, sidewalks, bicycles paths, etc. Thus, there 

7 exists a need to continue to collect information for a geographic database. 

8 One type of information that is useful to include in a geographic database is road 

9 grade information. One way of measuring road grade is with a vehicle having a 

10 mechanical road grade sensor traveling the road. One shortcoming of the mechanical 

1 1 road grade sensors is that they may provide inaccurate road grade measurements when 

12 traveling over bumps or a rough road. Additionally, road grade may be derived with 

13 altitude at points along the road determined using surveying equipment. However, 

14 surveying equipment may require that the equipment be mounted in a stationary position 

15 and therefore is not useable in a moving vehicle. 

16 Accordingly, it would be beneficial to collect and compute road grade information 

17 more inexpensively, efficiently and accurately. 
18 

19 SUMMARY OF THE INVENTION 

20 To address these and other objectives, the present invention comprises a method 

21 for obtaining road grade data for a geographic database. Data indicating a plurality of 

22 locations along roads including data indicating altitude at the plurality of locations is 

23 provided. Road grade values for a plurality of subsections of the road are computed; each 

24 of the subsections is a portion of the road between two of the locations. A plurality of 

25 temporary change points are identified as joining two adjacent subsections having road 

26 grade values that differ by a value greater than a predetermined amount. A road grade 

27 value for each part of the road between adjacent temporary change points is computed. 

28 Consecutive parts of the road having road grade values that falls within a predetermined 

29 range are grouped into at least one constant road grade section. The geographic database 

30 stores data that define the constant road grade section. 
31 
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1 BRIEF DESCRIPTION OF THE DRAWINGS 

2 Figure 1 shows a map of a geographic region. 

3 Figure 2 is a block diagram of a geographic database that represents the 

4 geographic region of Figure 1 . 

5 Figure 3 is a block diagram showing a road segment data record contained in the 

6 geographic database of Figure 2. 

7 Figure 4 is a block diagram illustrating a vehicle equipped for data collection. 

8 Figure 5 is a flow chart of the operations of a data collection facility. 

9 Figure 6 is a flow chart of the operations of the data collection facility for 

10 analyzing data to obtain road grade data. 

1 1 Figure 7 is a graph of altitude data points. 

12 Figure 8 is a graph of smooth altitude profile. 

13 Figure 9 is a graph of smooth altitude profile with grade change points. 

14 Figure 10 is a graph of smooth altitude profile with sections of constant road 

15 grade. 

16 Figure 1 1 is a graph that illustrates the Line Hough Transform. 

17 Figures 12A-12D are more graphs that illustrate the Line Hough Transform. 

18 Figure 13 is a flow chart of the operations of the data collection facility for using 

19 the Hough Line Transform to obtain road grade data. 

20 Figure 14 is a graph of altitude data analyzed using the Line Hough Transform. 
21 

22 DETAILED DESCRIPTION OF THE PRESENTLY PREFERRED EMBODIMENTS 

23 

24 I. GEOGRAPHIC DATABASE 

25 Figure 1 shows a map 10 of a geographic region 12. The geographic region 12 

26 may correspond to a metropolitan or rural area, a state, a country, or combinations 

27 thereof, or any other area of comparable size. Located in the geographic region 12 are 

28 physical geographic features, such as roads, points of interest (including businesses, 

29 municipal facilities, etc.), lakes, rivers, railroads, municipalities, etc. 

30 Figure 1 also includes an enlarged map 14 of a portion 16 of the geographic 

31 region 12. The enlarged map 14 illustrates part of the road network 18 in the geographic 

32 region 12. The road network 18 includes, among other things, roads and intersections 
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1 located in the geographic region 12. As shown in the portion 16, each road in the 

2 geographic region 12 is composed of one or more road segments 20. A road segment 20 

3 represents a portion of the road. Each road segment 20 is shown to have associated with 

4 it two nodes 22; one node represents the point at one end of the road segment and the 

5 other node represents the point at the other end of the road segment. The node at either 

6 end of a road segment may correspond to a location at which the road meets another road, 

7 i.e., an intersection, or where the road dead-ends. 

8 Referring to Figure 2, a geographic database 30 contains data 32 that represents 

9 some of the physical geographic features in the geographic region (12 in Figure 1). The 

10 data 32 contained in the geographic database 30 includes data that represent the road 

1 1 network 18. In the embodiment of Figure 2, the geographic database 30 that represents 

12 the geographic region 12 contains at least one database record 34 (also referred to as 

13 "entity" or "entry") for each road segment 20 in the geographic region 12 in Figure 1. 

14 The geographic database 30 that represents the geographic region 12 also includes a 

15 database record 36 (or "entity" or "entry") for each node 22 in the geographic region 12. 

16 (The terms "nodes" and "segments" represent only one terminology for describing these 

17 physical geographic features and other terminology for describing these features is 

18 intended to be encompassed within the scope of these concepts.) 

19 The geographic database 30 may also include other kinds of data 40. The other 

20 kinds of data 40 may represent other kinds of geographic features or anything else. The 

21 other kinds of data may include point of interest data. For example, the point of interest 

22 data may include point of interest records comprising a type (e.g., the type of point of 

23 interest, such as restaurant, hotel, city hall, police station, historical marker, ATM, golf 

24 course, etc.), location of the point of interest, a phone number, hours of operation, etc. 

25 The geographic database 30 also includes indexes 42. The indexes 42 may include 

26 various types of indexes that relate the different types of data to each other or that relate 

27 to other aspects of the data contained in the geographic database 30. 

28 Figure 3 shows some of the components of a road segment data record 34 

29 contained in the geographic database 30. The road segment record 34 includes a segment 

30 ID 34(1) by which the data record can be identified in the geographic database 30. Each 

31 road segment data record 34 has associated with it information (such as "attributes", 
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1 "fields", etc.) that describes features of the represented road segment. The road segment 

2 data record 34 may include data 34(2) that indicate the restrictions, if any, on the 

3 direction of vehicular travel permitted on the represented road segment. The road 

4 segment data record 34 includes data 34(3) that indicate a speed limit or speed category 

5 (i.e., the maximum permitted vehicular speed of travel) on the represented road segment. 

6 The road segment data record 34 may also include data 34(4) indicating whether the 

7 represented road segment is part of a controlled access road (such as an expressway), a 

8 ramp to a controlled access road, a bridge, a tunnel, a toll road, a ferry, and so on. 

9 The road segment record 34 may also include road grade data 34(5) that indicate 

10 the grade of the road segment. In one embodiment, the road grade data 34(5) includes 

1 1 road grade change points and a corresponding percentage of grade change. Additionally, 

12 the road grade data 34(5) may include the corresponding percentage of grade change for 

13 both directions of a bi-directional road segment. The location of the road grade change 

14 point is represented as a position along the road segment, such as thirty feet from the end 

15 or node of the road segment. For example, the road segment may have an initial road 

16 grade associated with its beginning node. The road grade change point indicates the 

17 position on the road segment wherein the road grade or slope changes, and percentage of 

18 grade change indicates a percentage increase or decrease of the grade or slope. Each road 

19 segment may have several grade change points depending on the geometry of the road 

20 segment. In another embodiment, the road grade data 34(5) includes the road grade 

21 change points and an actual road grade value for the portion of the road segment after the 

22 road grade change point until the next road grade change point or end node. In a further 

23 embodiment, the road grade data 34(5) includes elevation data at the road grade change 

24 points and nodes. 

25 The road segment data record 34 also includes data 34(6) providing the 

26 geographic coordinates (e.g., the latitude and longitude) of the endpoints of the 

27 represented road segment. In one embodiment, the data 34(6) are references to the node 

28 data records 36 that represent the nodes corresponding to the endpoints of the represented 

29 road segment. The road segment data record 34 may also include or be associated with 

30 other data 34(7) that refer to various other attributes of the represented road segment. 

31 The various attributes associated with a road segment may be included in a single road 
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1 segment record, or may be included in more than one type of record which are cross- 

2 referenced to each other. For example, the road segment data record 34 may include data 

3 identifying what turn restrictions exist at each of the nodes which correspond to 

4 intersections at the ends of the road portion represented by the road segment, the name or 

5 names by which the represented road segment is known, the street address ranges along 

6 the represented road segment, and so on. 

7 Each of the node data records 36 may have associated information (such as 

8 "attributes", "fields", etc.) that allows identification of the road segment(s) that connect 

9 to it and/or its geographic position (e.g., its latitude and longitude coordinates). For the 

10 embodiment shown in Figure 3, the node data records 36(1) and 36(2) include the latitude 

1 1 and longitude coordinates 36(1)(1) and 36(2)(1) for their node. The node data records 

12 36(1) and 36(2) also include other data 36(1)(2) and 36(2)(2) for their node. 

13 The data records 34 in the geographic database 30 that represent roads may not 

14 necessarily include all the same types of data attributes. One reason for this is that roads 

15 do not all have the same properties. For example, some roads have a highway 

16 designation (e.g., "Wisconsin State Highway 120") whereas other roads do not. Another 

17 reason why data records in the geographic database 30 that represent roads may not have 

18 the same data attributes is that some of the properties of a road may not have been 

19 collected or confirmed. 

20 Road grade data may be collected using the data collection system as described 

21 below. In one embodiment, road grade data are obtained for all the roads represented in 

22 the geographic database. In another embodiment, road grade data are included for only 

23 some of the roads represented in the geographic database. According to this latter 

24 embodiment, some of the roads are represented by data records that do not include road 

25 grade data. 



26 

27 H. SYSTEM FOR OBTAINING ROAD GRADE DATA 

28 A. Data Collection System 

29 Figure 4 shows a vehicle 50 equipped with a data collection unit 52. The vehicle 



30 50 may be a data collection vehicle operated by a geographic database developer. The 

31 data collection vehicle is operated specifically for data collection purposes. The 



6 



N0155US 



1 geographic database developer determines the routes traveled by the vehicle in order to 

2 collect specific data. Alternatively, the vehicle 50 is operated as a probe vehicle. A 

3 probe vehicle is a vehicle that collects road-related data while it is being used for 

4 purposes unrelated to the collection of road-related data. For example, a probe vehicle is 

5 operated for ordinary, everyday purposes, such as commuting, leisure or business. A 

6 member of the public may operate the probe vehicle, or alternatively, a commercial 

7 enterprise or government entity may operate the probe vehicle. While the probe vehicle 

8 is being used for a purpose unrelated to the collection of road-related data, the data 

9 collection unit 52 in the vehicle 50 operates to collect road-related data. 

10 The data collection unit 52 is a combination of hardware and software 

11 components. The data collection unit 52 includes a positioning system 54. The 

12 positioning system 54 includes a GPS 56. The positioning system 54 may also include 

13 inertial sensors 58. The positioning system 54 enables the position of the data collection 

14 unit 52 (and thus the vehicle 50 in which it is located) to be determined. The positioning 

15 system 54 provides a coordinate triple of latitude, longitude, and altitude. In one 

16 embodiment, the GPS 56 has an accuracy of approximately 3 meters for latitude and 

17 longitude; however, the GPS 56 may have any other accuracy, such as accuracy ranging 

18 from 2 centimeters to 20 meters for latitude and longitude. The uncertainty in GPS 

19 generated altitude is significantly worse than for latitude and longitude. The data 

20 collection unit 52 also includes other sensors 60, such as a speed sensor indicating the 

21 current speed of the vehicle 50, an odometer indicating distance traveled, a sensor that 

22 indicates heading or bearing of the vehicle or any other sensor. In another embodiment, 

23 the positioning system 54 may include one sensor to provide latitude and longitude 

24 position and another separate sensor to provide altitude. 

25 The data collection unit 52 also includes the necessary hardware and software 

26 (processor and programming 62) to receive data from the positioning system 54 and other 

27 sensors 60. The data collection unit 52 includes the appropriate interfaces to allow the 

28 processor and programming 62 to receive data from the positioning system 54 and the 

29 other sensors 60. The processor and programming 62 in the data collection unit 52 are 

30 suitable for selecting or receiving the position data from the positioning system 54 at 

31 appropriate time periods, such as every second. Additionally, the processor and 
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1 programming 62 select or receive data from the other sensors 60 at corresponding times. 

2 The processor and programming 62 relate the data received from the positioning system 

3 54 and the other sensors 60 to each other. Furthermore, the processor and programming 

4 62 may associate the data received from the positioning system 54 and the other sensors 

5 60 with a time stamp or any other information. 

6 The data collection unit 52 includes a data storage device 64. The processor and 

7 programming 62 in the data collection unit 52 provide for storing some or all the data 

8 from the positioning system 54 and the other sensors 60 on the data storage device 64. 

9 The data storage device 64 is a non-volatile data storage unit, such as a hard drive or a 

10 memory card. 

11 In one embodiment, the data collection unit 52 may receive data from a vehicle 

12 navigation system or may be part of the navigation system. There are various kinds of 

13 navigation systems installed in vehicles. Some navigation systems use vehicle data in 

14 connection with providing certain navigation-related functions. As an example, some 

15 vehicle navigation systems use the vehicle speed or wheel speed in conjunction with 

16 other devices, such as a gyroscope or GPS, to help determine the position of the vehicle. 

17 The data collection unit 52 may receive position and speed data from the navigation 

18 system. Alternatively, the navigation system may perform the functions of the data 

19 collection unit 52 and may include an additional processor and programming to collect 

20 data. 

21 The data collection unit 52 may include a communications system 66. The 

22 communications system 66 provides for sending some or all the data from the positioning 

23 system 54 and the other sensors 60 to a remotely located data collection facility 68. In 

24 the embodiment of Figure 4, the communications system 66 provides for wireless 

25 transmission of the data from the data collection unit 52 to the remotely located data 

26 collection facility 68. In the embodiment of Figure 4, the data sent from the data 

27 collection unit 52 are stored temporarily on the data storage device 64 before being 

28 transmitted to the data collection facility 68. Various processes may be performed on the 

29 data before they are sent to the remotely located data collection facility 68. For example, 

30 the data may be compressed, filtered, normalized, etc. These processes may be 
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1 performed to reduce the amount of data that need to be sent from the data collection unit 

2 52 to the data collection facility 68. 

3 In one alternative, the data collected by the data collection unit 52 are sent every 

4 several seconds or minutes to the data collection facility 68 after being temporarily stored 

5 on the data storage device 64 in the data collection unit 52. In another alternative 

6 embodiment, the data collected by the data collection unit 52 are sent directly to the data 

7 collection facility 68 without being temporarily stored on a data storage device in the data 

8 collection unit 52. In another alternative embodiment, the data collected by the data 

9 collection unit 52 are stored in the vehicle for a relatively long period of time, e.g., 

10 several days or weeks. The data are then sent to the data collection facility 68. The data 

1 1 may be sent wirelessly via a communications system or alternatively, the data may be 

12 sent by other means. For example, the data may be transmitted over land-based 

13 telephone lines or the Internet. In another embodiment, the data storage medium upon 

14 which the data are stored is physically sent to the data collection facility 68 (e.g., by 

15 mail). At the data collection facility 68, further processing of the data takes place as 

16 described below. 

17 B. Data Collection Facility 

18 According to one embodiment, the data collection facility 68 obtains position 

19 data, including altitude data, computes road grade information and updates the 

20 geographic database following the steps shown in Figure 5. The data collection facility 

21 68 includes necessary hardware and software (processor and programming) to carry out 

22 the steps of Figure 5. The data collection facility 68 receives the data from the vehicle 50 

23 at step 70. The data collection facility 68 obtains the data from the data collection 

24 vehicle that is operated by the geographic database developer. Alternatively, the data 

25 collection facility 68 may obtain the data from numerous probe vehicles. The data 

26 collection facility 68 may obtain the data from the vehicle(s) 50 by wireless data 

27 transmission or by other means (e.g., sending a diskette or via modem). 

28 The data collection facility 68 processes the data received from the vehicle 50 and 

29 computes road grade data at step 72. The methods for computing road grade data will be 

30 described in detail below in conjunction with Figures 6 through 12. 
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1 Referring to Figure 5, after the computing of step 72, the road grade data are 

2 stored in a master copy 78 of the geographic database at step 74. The road grade data 

3 may be stored as an attribute to the road segment data record 34. The road grade data 

4 that are stored in the master copy 78 of the geographic database may be used to update 

5 existing data or to add new data. For example, the master copy 78 of the database may 

6 already include road grade data for a particular represented road segment. The new road 

7 grade data obtained by using the process described in Figure 5 can be used to update the 

8 existing data, e.g M confirm the existing data or make the existing data more accurate. 

9 Alternatively, the master copy 78 of the geographic database may not include road grade 

10 data for a particular road segment. If new road grade data are obtained for a road 

1 1 segment, the new road grade data can be added as a new attribute of the data record. 

12 In one embodiment, prior to updating existing data or adding new data to the 

13 master copy of the geographic database with the road grade data, the existence and 

14 position of some the identified road grade information is confirmed by direct observation 

15 of field personnel. 

16 The geographic database with new or improved road grade data and other data 

17 can be used to make derived database products at step 76. The derived database products 

18 may include only portions of all the data in the master version 78 of the database. For 

19 example, the derived database products may include data that relate to only one or more 

20 specific regions. The derived database products may be used on various kinds of 

21 computing platforms. For example, the derived database products may be used in 

22 navigation systems (such as in- vehicle navigation systems and hand-held portable 

23 navigation systems), personal computers (including desktop and notebook computers), 

24 and other kinds of devices (such as PalmPilot®-type devices, pagers, telephones, 

25 personal digital assistants, and so on). Derived database products may also be used on 

26 networked computing platforms and environments, including the Internet. The derived 

27 database products may be in a different format than the format in which the master copy 

28 of the database is maintained. The derived database products may be in a format that 

29 facilitates the uses of the derived products in the platforms in which they are installed. 

30 The derived database products may also be stored in a compressed format on the media 

31 on which they are located. 
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1 In an alternative embodiment, the data collection unit 52 of Figure 4 may perform 

2 some of the steps of the data collection facility 68. In this embodiment, rather than 

3 sending the position data, including altitude data, to the data collection facility 68, the 

4 processor and programming 62 of the data collection unit 52 performs the analysis on the 

5 data to obtain the road grade information. The data collection unit 52 then stores the road 

6 grade data in a storage medium. 
7 

8 HI. METHOD FOR COMPUTING ROAD GRADE DATA 

9 Figure 6 illustrates one embodiment of the steps performed by the data collection 

10 facility 68 to compute road grade data of step 76 in Figure 5. Because the positioning 

11 system 54 used to collect the altitude and latitude/longitude data may have accuracy 

12 limitations, the road grade determination may be complicated by an uneven distribution 

13 of data points (coordinate triple of latitude, longitude and altitude) and the presence of 

14 errors in the altitude and latitude/longitude data. The uneven distribution of data and 

15 errors in the position data may suggest that the road grade zigzags between inclines and 

16 declines instead of the gradual road grade changes present in typical road geometry. To 

17 reduce the effect of the errors in the altitude and latitude/longitude data and the uneven 

18 distribution of data points, the data collection facility 68 filters the data. At step 80, the 

19 data collection facility 68 filters the data to remove non-uniform and closely spaced data 

20 points. When the vehicle 50 is stopped or moving at slow speeds provides a large 

21 number of closely spaced data points. Due to the accuracy of the positioning system, the 

22 closely spaced data points may provide an inaccurate representation of a true altitude 

23 profile of the road. 

24 Figure 7 illustrates an altitude profile of altitude data versus horizontal distance 

25 along the road as collected by the vehicle 50 traveling along the road. Horizontal 

26 distance is the accumulated distance between latitude/longitude positions of consecutive 

27 data points along the road segment. Referring to Figure 7, box 92 illustrates a cluster of 

28 closely spaced data points 94. In one embodiment, the data collection facility 68 

29 eliminates the cluster of closely spaced data points 94. In another embodiment, the data 

30 collection facility 68 replaces cluster of closely spaced points with a linear representation, 

31 such as the dashed line 96 in Figure 7. 
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1 After eliminating the closely spaced data points, the data collection facility 68 

2 derives a smooth altitude profile of the road segment at step 82 of Figure 6. Even after 

3 eliminating closely spaced data points, the data points may provide an uneven altitude 

4 profile 98 as shown in Figure 8. Since roads by design are smooth surfaces with gradual 

5 altitude changes, the data collection facility 68 derives the smooth altitude profile to 

6 reduce such unrealistic profile 98. In one embodiment, the data collection facility 68 uses 

7 a least squares polynomial algorithm to derive the smooth altitude profile; however, other 

8 smoothing or curve fitting algorithms may also be used. 

9 In one embodiment, the data collection facility 68 derives the smooth altitude 

10 profile using a least squares polynomial with (2n+l) successive altitude data points. In 

11 the least squares polynomial implementation, a fourth order polynomial uses successive 

12 altitude data points to determine the coefficients of the polynomial (Equation 1). 

13 
14 
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17 As shown in the Equation 1, the independent variable is the horizontal distance (d) from 

18 the first to each of the following latitude/longitude points and the dependent variable is 

19 their corresponding altitude difference (A/i). In Equation 1, the coefficients of the 

20 polynomial are represented by variable a. Equation 2 represents Equation 1 in matrix 

21 notation: 
22 

23 y = Ax (2) 

24 

25 Where, x = [a 0 a x a 2 a 3 a 4 f (T is symbol for transposed matrix). 

26 Using least squares principle for the (2n+l) equations, the five parameters 

27 (polynomial coefficients) are optimally determined by matrix methods (linear algebra) 

28 such as the pseudo inverse or from the normal equation (Equation 3). 
29 
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1 x =[a t a)~ 1 A T y (3) 

2 

3 Once these parameters are determined, new value for the altitude of the middle (n+l*) 

4 point is obtained from Equation 4. Altitude difference (original altitude - new altitude) 

5 for each point is also derived (Equation 5). 
6 

7 h n = A, +a 0 + a x d n + a x d 2 n + a x d\ + a x d A n (4) 

8 

9 Ah n =h n - h n (5) 

10 

11 Once new altitude and altitude difference for a point are obtained, the (2n+l) window is 

12 shifted one point forward and the process is repeated until the end of the profile. Once 

13 this process is complete, variance (a 2 ) of the altitude difference is calculated using 

14 Equations 6 and 7. 
15 

16 M = ^—^ (6) 



17 

18 a\ = ^ V 1 (7) 

n-1 

19 
20 

21 Large variance values obtained from Equation 7 indicate points having larger 

22 altitude errors. If the variance value is larger than a chosen threshold, such as one or two 

23 standard deviations (standard deviation is square root of variance), an outlier point is 

24 identified. Figure 8 illustrates some outlier points 100. 

25 The least squares polynomial process is then repeated to re-estimate the altitudes 

26 excluding the identified outlier data points 100. The removal of outlier data points 100 

27 ensures that incorrect altitude data will not influence the final smooth altitude profile. In 

28 one embodiment, the smoothing process is also repeated in reverse order to get a second 

29 value of altitude for the each point. A mean altitude value calculated from the forward 

30 and reverse calculated altitude values further reduces the least square polynomial 

31 estimation errors. In one embodiment, the above process is run iteratively until a 
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1 significantly smooth altitude profile is obtained. Figure 8 illustrates the smooth altitude 

2 profile 102 provided by the least squares polynomial process. 

3 In one embodiment, after deriving the smooth altitude profile, the data collection 

4 facility 68 calculates road grade values between successive data points (sub-segments of 

5 the road segment) at step 84 of Figure 8. The road grade values for each of the sub- 

6 segments is calculated as the change in altitude divided by the change in horizontal 

7 distance between the two successive data points. Road grade is typically expressed in 

8 percentage, so the previously calculated change in altitude divided by change in 

9 horizontal distance is converted into a percentage. The change in altitude is calculated 

10 using the smoothed altitude data associated with each data point, and the change in 

1 1 horizontal distance is calculated using the latitude/longitude data associated with each 

12 data point. 

13 Once the road grade values for each sub-segment are calculated, the data 

14 collection facility 68 groups together a number of adjoining sub-segments with grades 

15 values within a selected range and categorizes these adjoining sub-segments into a 

16 section of constant road grade. According to one embodiment, the data collection facility 

17 68 derives the sections of constant road grade following steps 86, 88 and 90 of Figure 6. 

18 At step 86, the data collection facility 68 identifies zero-crossing points and pseudo zero- 

19 crossing points in the smooth altitude profile 102. Zero-crossing points correspond with 

20 data points that join successive sub-segments whose road grade values change from 

21 positive to negative or vice- versa. Zero-crossing points are critical locations of peaks and 

22 valleys on the road segment corresponding to road grade changes. Moving along the 

23 road segment and comparing the road grade value of one of the sub-segments to the road 

24 grade value of the next sub-segment identifies the location of zero-crossing points as the 

25 data point corresponding with the grade values changing from negative to positive or 

26 vise-versa. Figure 9 illustrates zero-crossing points 104 on the smooth altitude profile 

27 102. 

28 Zero-crossing points may not include all locations on the road where road grade 

29 changes because significant road grade changes also occur between consecutive zero- 

30 crossing points. For example, the road may change from a highly steep road grade to a 

31 significantly less steep road grade. In one embodiment, to find locations corresponding 
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to road grade changes between the identified zero-crossing points, the data collection 
facility 68 implements a rotational transformation. The rotational transformation rotates 
the original altitude profile between two consecutive zero-crossing points by a chosen 
angle. The rotation generates a pseudo profile with the zero-crossing points turned to 
non-zero-crossing points and any points where the slope changes significantly are false 
zero-crossing points. These false zero-crossing points obtained with the rotation are 
called pseudo zero-crossing points. 

The process for detecting pseudo zero-crossing points first calculates the altitude 
difference (Afi) and the distance difference (D) to all the points from the first zero- 
crossing point to the second zero-crossing points with Equations 8 and 9 wherein (Ek - 
Ei) represent displacement in easting and (N* - NO represent displacement in northing 
from point 1 to point k. A slope angle (0) for the stretch between the two zero-crossing 
points is then obtained using Equation 10. A fraction of this angle (multiplied by the 
fractional constant k), for example k equals one half, is used to determine a rotation 
matrix R with Equation 11. 
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1 Using the rotation matrix with the calculated altitude and distance differences, pseudo- 

2 altitudes for all the points between the two zero-crossing points are calculated with 

3 Equation 12. Using the pseudo altitude profile provided with pseudo-altitude calculation, 

4 pseudo zero-crossing points are identified as data points that join successive sub- 

5 segments whose road grade values change from positive to negative or vice-versa. Figure 

6 9 illustrates pseudo zero-crossing points 106 on the smooth altitude profile 102. 

7 In an alternative embodiment, rather than identifying the zero-crossing points and 

8 pseudo zero-crossing points to identify locations on the road where road grade changes 

9 significantly in the manner described above, the data collection facility 68 uses another 

10 method to identify any data points that join successive sub-segments whose road grade 

1 1 values differ by more than a predetermined amount as a location on the road where the 

12 road grade significantly changes. 

13 The identified zero-crossing points and pseudo zero-crossing points provide 

14 points identifying locations on the road where the road grade significantly changes. In 

15 one embodiment, the data collection facility 68 performs additional processing on the 

16 identified grade change points to remove unnecessary grade change points due to residual 

17 inaccuracy in the altitude and the latitude/longitude data. At step 88 of Figure 6, the data 

18 collection facility 68 calculates the road grade value between grade change points and 

19 categorizes the sections between grade crossing points into grade bins. The grade bins 

20 represent a range of road grade values. Table 1 illustrates one example of the grade bins 

21 and their respective road grade value ranges; other ranges of grade values are also 

22 possible. 

23 Table 1 



Grade Bin 


Grade Range 


Grade Bin 


Grade Range 


1 


grade < -15 


10 


1 <= grade < 3 


2 


-15 <= grade < -13 


11 


3 <= grade < 5 


3 


-13 <= grade < -11 


12 


5 <= grade < 7 


4 


-11 <= grade < -9 


13 


7 <= grade < 9 


5 


-9 <= grade < -7 


14 


9 <= grade < 1 1 


6 


-7 <= grade < -5 


15 


11 <= grade < 13 


7 


-5 <= grade < -3 


16 


13 <= grade < 15 


8 


-3 <= grade < -1 


17 


grade >= 15 


9 


-1 <= grade < 1 
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1 Moving along the road segment, the data collection facility 68 calculates the road 

2 grade value from the first grade change point to the second grade change point and the 

3 road grade value from the first grade change point to the third grade change point. If the 

4 two calculated road grade values fall into the same grade bin, the section from the first to 

5 the third grade change point forms one section of constant road grade, and the second 

6 grade change point is eliminated. The process continues forward to the fourth grade 

7 change point with the road grade value from the first to the fourth grade change point 

8 being compared to the grade value from the first to the third grade change point. If the 

9 two calculated road grade values fall within the same grade bin, the process eliminates 

10 the third grade change point and continues forward to the next grade change point. 

1 1 However, if the two road grade values do not belong to the same grade bin, the process 

12 retains the third grade change point, and the section of constant road grade extends from 

13 the first to the third road grade change point. The above procedure repeats from the third 

14 grade change point forward to the next grade change point to identify the next section of 

15 constant road grade. The process proceeds until reaching the forward most grade change 

16 point on the road segment. In addition to removing unnecessary grade change points, the 

17 data collection facility 68 identifies and adds locations for additional grade change points 

18 when the road grade value of one section differs greatly from the road grade value of the 

19 adjacent section, such as one falls within grade bin 3 and the other falls within grade bin 

20 8. 

21 Residual inaccuracies in the collected data may lead to significant road grade 

22 value changes for short sections of the road segment than actually present on the road 

23 segment. Because roads are designed with certain constraints to allow safe driving, large 

24 road grade value changes unlikely occur within short sections of the road segment. 

25 Therefore, the data collection facility 68 filters the sections of constant road grade using a 

26 minimum grade distance for road grade changes and merges short identified sections of 

27 constant road grade at step 90 of Figure 6. The minimum grade distance may be based on 

28 road design criteria, truck length, type of terrain, etc. The data collection facility 68 

29 compares the road length of identified sections of constant road grade to the minimum 

30 grade distance. If any of the identified sections of constant road grade are shorter than 

31 the minimum grade distance, the short section is merged with either the preceding or 
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1 following adjoining section having the most compatible road grade value. When all of 

2 the short sections are merged, the remaining grade change points closely model the actual 

3 road geometry. 

4 Also at step 90 of Figure 6, the data collection facility 68 identifies any abnormal 

5 road grade values between the remaining grade change points. Typically, road geometry 

6 design allows a restricted range of road grade values for different types of roads and/or 

7 terrain. If any of the road grade values between the grade change points do not fall 

8 within the range of values based on road geometry design, such sections are merged with 

9 the adjacent section having the closest road grade value and the grade change point 

10 between the sections is eliminated. This process is repeated until all grades are within 

11 admissible range. In one embodiment, the admissible range for road grade values is 

12 between negative and positive ten percent. In another embodiment, the admissible range 

13 for road grade values may vary depending on the landscape and information regarding 

14 road construction. In an alternative embodiment, the abnormal grade values can be 

15 reported for human interpretation and/or field inspection. 

16 Figure 10 illustrates the identified grade change points 108 and sections of 

17 constant road grade 1 10 for a road segment as determined by the above method. The 

18 location of the grade change points 108 and the road grade values corresponding to the 

19 sections of constant road grade 1 10 may be included as the road grade data 34(5) of the 

20 road segment data record 34 as shown in Figure 3. Additionally, the direction of the road 

21 grade may be determined for inclusion with the road grade data. The direction of road 

22 grade for the sections of constant road grade may be determined by comparing the 

23 altitude value associated with the road grade change points that define the section of 

24 constant road grade. For example, if one grade change point has a corresponding altitude 

25 value greater than the altitude value of the next grade change point indicates an incline, if 

26 one grade change point has a corresponding altitude value less than the altitude value of 

27 the next grade change point indicates a decline, and if one grade change point has a 

28 corresponding altitude value equal to the altitude value of the next grade change point 

29 indicates a flat section. Additionally, once road grade values have been determined for 

30 each of the sections of constant road grades, percentage change of grade from one section 

31 of constant road grade to an adjacent section of constant road grade may be calculated. 
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1 Although the above description of computing the road grade data has been 

2 utilizing altitude data collected with the positioning system 54, the above method may 

3 also be applied to other collected data. For example, the method may determine road 

4 grade data using altitude data collected with a barometer as described in the co-pending 

5 application entitled "METHOD AND SYSTEM FOR OBTAINING ROAD GRADE 

6 DATA" filed on the same date herewith, Attorney Docket No. N0154US, the entire 

7 disclosure of which is incorporated by reference herein. Additionally, the method may 

8 determine road grade data using height data from digital elevation model, pitch from a 

9 gyroscope, height difference from a double GPS antenna, and forward and vertical (delta) 

10 velocity from accelerometers or GPS. 
n 

12 IV. ALTERNATIVES 



13 A. Hough Transform 

14 In an alternative embodiment, the data collection facility 68 implements a Hough 

15 Transform on the altitude and latitude/longitude data points to identify the road grade 

16 data. The Hough Transform is a known technique for detecting straight lines or simple 

17 curves (such as circular arcs) within an image. (It can also be extended to other shapes, 

18 to some extent). The Hough Transform is described in U.S. Pat. No. 3,069,654. 

19 In the case of line detection, the Hough Transform will find all parameter pairs 

20 (ah,c) that describe each line segment y = mx + c occurring in the image, where m is 

21 slope and c is the y -intercept. Note that for any point (x,y) in the x-y plane, an 

22 infinite number of lines can be drawn through that point, each with a different slope and 

23 y -intercept as shown in Figure 1 1. 

24 The equation of any one of these lines is 

25 y = mx + c (1) 

26 However, since m and c can take on a value of ±«> , it is more convenient to 

27 write the equation of the line in parametric form as follows: 

28 p = xcos0 + ysin0 (2) 
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1 where 9 is the angle made with the x -axis by the normal through the origin to the line 

2 and p is the length of the normal, as shown in Figures 12A, 12B, 12C and 12D. The 



3 



10 



relationship between {p,0) and (m,c) is given below: 



0 = ^-tan l (m) (3) 



Vm 2 +1 



6 Each parameter pair (/?,#) defines a unique line in the x - y plane. Each point in the jc - 

7 y plane maps to a sinusoid in the/? - 0 plane, i.e., if we plot in the p - 0 plane, each 

8 pair (p,0) which defines a line through the point (x, y) , the resulting curve is a sinusoid. 

9 The Line Hough Transform of a point in the x - y plane is thus a sinusoid in the 
p - 0 plane. Points in the p - 0 plane where sinusoids intersect correspond to the 

11 existence of a single line in the x - y plane. In practice, the intersection points are 

12 detected using an accumulator, which divides the (x - a) 2 + (y - b) 2 = r 2 plane into a 

13 grid. Each point in the x - y plane produces a vote in a cell in the accumulator. Cells 

14 with a large number of votes correspond to an intersection point, which means that a line 

15 segment is present in the image in the x - y plane. 

16 Figure 13 illustrates the steps performed by the data collection facility 68 to 

17 implement the Hough Transform on the altitude and latitude/longitude data to identify 

18 road grade data. At step 120, the data collection facility 68 computes the Hough Line 

19 Transform with the horizontal distance along the road as the x-axis variable and the 

20 corresponding altitude data as the y-axis variable. Horizontal distance is the accumulated 

21 distance between consecutive latitude/longitude data points along the road. 

22 As discussed above, the Hough Line Transform identifies lines in the data. At 

23 step 122, the data collection facility 68 identifies lines from the Hough Line Transform as 

24 sections of constant road grade along the road. Figure 14 illustrates a graph of altitude 

25 versus horizontal distance along a road. The sections of constant road grade correspond 

26 to the lines 128 identified with the Hough Line Transform. 
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1 At step 124, the data collection facility 68 identifies road grade change points 130 

2 at intersections of adjacent identified lines from the Hough Line Transform. The grade 

3 change points 130 in Figure 14 illustrate locations at which the road grade changes from 

4 one section of constant road grade 128 to another section of constant road grade 128. In 

5 one embodiment, the location of the grade change points correspond to the latitude and 

6 longitude position associated with the intersection point of the adjacent lines. 

7 At step 126, the data collection facility 68 identifies the road grade values of the 

8 sections of constant road grade as the slopes of the identified lines from the Hough Line 

9 Transform. The Hough Line Transform provides the slope of each line representing a 

10 section of constant road grade. The road grade value for each section of constant road 

1 1 grade is the slope value of its corresponding line. Since road grade is typically expressed 

12 in percentage, the slope value is converted into a percentage. The direction of the slope 

13 of the lines provides the direction of the road grade. For example, a positive slope 

14 indicates an incline, a slope of zero represents a flat section, and a negative slope 

15 indicates a decline. In another embodiment, the direction of road grade may be 

16 determined by comparing the altitude value associated with road grade change points. 

17 For example, if one grade change point has a corresponding altitude value greater than 

18 the altitude value of the next grade change point indicates an incline. 
19 

20 B. Applications for the Road Grade Data 

21 The road grade data within the geographic database have several applications. 

22 First, various driver assistance systems for vehicles that help make travel more efficient 

23 and convenient may use the road grade data. For example, one driver assistance system 

24 that may use the road grade data is an automatic transmission control system. The 

25 automatic transmission control system adjusts the operation of the vehicle transmission to 

26 match the grade of the upcoming road. Another driver assistance system that may use the 

27 road grade data is an adaptive cruise control system. 
28 

29 It is intended that the foregoing detailed description be regarded as illustrative 

30 rather than limiting and that it is understood that the following claims including all 

31 equivalents are intended to define the scope of the invention. 
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